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Abstract 

We study a boundary perturbation problem for a one dimensional Schrodinger equation 
in which the potential has a regular singularity near the perturbed end point. We give the 
<**] ■ asymptotic behaviour of the eigenvalues under the perturbation. This problem arose out 

"£i ' of the author's studies of singular elliptic operators in higher dimensions and we illustrate 

this point with an example. The class of potentials to which this method applies is larger 
than that covered by standard results which assume uniform ellipticity of the operator or a 
perturbative term which is analytic in the perturbation parameter. 
£\ AMS Subject Classifications:47E05, 34L20, 34L40, 34L99, 65L15. 

^ ' Keywords:one dimensional Schrodinger operator, spectral convergence, asymptotic theory 

^\ , of eigenfunctions, domain perturbation. 
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^ ■ 1 Introduction 

O 

This paper forms part of the author's studies of the behaviour of the spectrum of the Dirichlet 
,S^ ■ Laplacian under perturbations of the domain. We let f2 be some bounded domain (either in R N 

or more generally in an incomplete Riemannian manifold) and then consider a slightly smaller 
region f2 e such that 

{x e £1 : dist(s, dty < e} C fi e C 

for e > small enough. Let H denote an elliptic operator acting in L 2 (f2) with Dirichlet boundary 
conditions and suppose it has eigenvalues A„. Now let \ n , e denote the eigenvalues of H restricted 
to fl e . The problem then is to study the asymptotics of the quantity A„ je — A„ as e — > + . 

The specific model we shall consider in this paper is a one-dimensional Schrodinger equation 

-f'(x) + V(x)f(x) = Xf(x) 
a < x < b 

subject to Dirichlet boundary conditions. The potential is assumed to have a regular singularity 
at the end point x = a and thus we are led to consider the boundary perturbation problem with 



the perturbed domain (a + e, b). More precisely we assume that 

V{x){x-df/c^> 1 

as x — > a + for some —1/4 < c < 3/4. We do not however assume that V(x)(x — a) 2 is analytic and 
so we cannot employ standard techniques to determine series expansions of the eigenfunctions as 
described in pj for example. 

The stability of the spectrum under perturbations of the boundary has a long history. In par- 
ticular we mention the texts of Kato, ||, Section 6.5], Courant-Hilbert, [|| p419], and Reed and 
Simon, || Chapter XII]. These may be used to show that A„ i£ are holomorphic in e near e = 
provided the operator is uniformly elliptic and Q £ may be obtained from Q by a sufficiently well- 
behaved transformation. These techniques work by changing the problem from one in which the 
domain is a function of some small parameter e and the operator is fixed to one in which the 
domain is fixed and the operator depends on the small parameter in some analytic way. 

If the Euclidean region is star-shaped then the eigenvalues of the Laplacian are known to satisfy 

A n , e — A„ = 0(e) as e — ► + . 

However, in pj the author considered a completely general type of boundary perturbation as- 
suming no regularity of the coefficients or indeed of the boundary. It merely assumed that some 
generalised Hardy inequality was valid. We refer to this paper for details but stress that the 
results can be applied when the operator is singular, including operators of the form 

Hf := -d(x) a A E (1) 

where d(x) is the distance to the boundary, < a < 1 and A^ is the ordinary Euclidean 
Laplacian in K 2 . To show the optimality of the results we considered the operator (pi) acting in a 
rotationally invariant domain and then used separation of variables to reduce the problem to the 
type of one-dimensional problem considered in this paper. This specific problem is dealt with in 



Example 4.2 



The main result, contained in Theorem 3.1, is that we have the asymptotic expansion 

A„ )£ = A„ + c n e p + o(e p ) 

as e — > + where p := 1\J c + 1/4. This result is important because it shows how to perform a 
boundary perturbation on a non-uniformly elliptic equation using only the standard theory of 
ODEs. In particular the standard results of Kato cannot possibly be applied. Moreover, Example 



4.2 and its connection with (hi) show that the result has implications for higher dimensional 



singular problems. 

There are many good introductions to the theory of ordinary differential equations and in par- 
ticular the classification of the boundary. We will not reproduce this theory here. The reader 
who is unfamiliar with the terminology or needs a refresher course could look for example at the 
comprehensive review of the subject in M. 

Finally to fix notation we remark that we write 

f(x) ~ g(x) as x — * a+ 



to mean 

f{x) , 

2 Set-up of the Problem 

Let X be the operator acting in L 2 (a,b) given initially by 

Lf := -f"(x) + V(x)f(x) for / e C 2 c (a, b). 
We assume the following: 

1. V is continuous throughout (a, 6). 

2. V is twice diffcrcntiable near x = a. 

3. V satisfies 

c 

^ ( x ) ~ 7 ^ as n a + . 

(x — a) z 

where —1/4 < c < 3/4. Moreover, this behaviour may be differentiated in the sense that 
\V'\ = 0((x - a)- 3 ) and \V"\ = 0{{x - a)" 4 ). 

4. L is bounded below and its spectrum consists solely of simple eigenvalues. 

Condition (1) will ensure that initial value problem is well-posed and has solutions in the classical 
sense. Conditions (2) and (3) specify the singularity at the end point is regular and will allow 
us to use Liouvillc- Green theory to determine the asymptotic behaviour of the eigenfunctions. 
The range of values of c is not arbitrary. If c < 1/4 then the corresponding operator ceases to be 
bounded below and if c > 3/4 then the spectral properties change (see below). 

The simplicity of the spectrum is guaranteed by taking separated boundary conditions. The fact 
that L is bounded below allows us to take the Fricdrichs extension and so deal with a self-adjoint 
operator. We comment that the Hardy Inequality immediately gives us the fact that L is bounded 
below provided V_, the negative part of the potential, satisfies 

for all x G (a, b) where d(x) = min{|x — a\, \x — b\} (condition (3) does not, of course, imply this 
inequality) . 



Lemma 2.1 Suppose 



c 

V(x) ~ 7 rr: as X - 

(x — a) z 



where c > —1/4. Moreover, suppose that this asymptotic behaviour can be differentiated. Then 
we have a basis of solutions comprising /+ and /_ where 

f±(x) ~ (x - a) 5± v c+ 3 as x — » a + . 
Moreover, if —1/4 < c < 3/4 i/ien 6oi/i /+ and /_ belong to L 2 (a,S) for some S > a. 

Proof The behaviour of the eigenfunctions is proved using Liouville Green theory. See |)[ Chapter 
6, Section 3] for details. The convergence or otherwise of / \f±(x)\ 2 dx is determined by the 
convergence or otherwise of J (x — a) 1 V c+ *dx for d near enough to a. 

The solution / + is unique up to constant multiples and is refered to as the recessive solution. 
The solution f% is non-unique and is called the dominant solution. 

In the case that both /+ and /_ are square integrable near the end point and vanish only finitely 
often we must give a boundary condition (this is known as the 'Limit Circle Non Oscillatory' case, 
denoted LCNO). The problem of specifying the correct boundary condition so that the equation 
does in fact correspond to taking the Friedrichs extension of the operator L has been solved. We 
briefly outline the result here; for more information see M. 

The maximal domain for L on (a, b), denoted by M{a, b) is defined by 

M(a, b) := {/ : (a, b) -» C : /, /' e AC{a, b) and 
f,Lf€L 2 (a,b)} 

where AC (a, b) denotes the space of absolutely continuous functions defined on (a, b). We define 
the Wronskian of u,v G M (a, S] by 

[w,w](a;) := — p(x)u'(x)v(x) + u(x)p(x)v' (x) 

and we make the definition 

[u,i>](a) :— lim [u, v](x). 

x — >a 

The correct LCNO Friedrichs boundary condition is then 

[f,u](a) = Q 

where u is a solution that satisfies (w 2 ) -1 €" L 2 (a, c) (u is called a principal solution at a and in 
the non-oscillatory case is just the recessive solution). 

Corollary 2.2 Let f + be the recessive solution near x = a. Then, in the case —1/4 < c < 3/4 
the Friedrichs boundary condition is given by 

[/,/+](«) = 

If the constant c is equal to —1/4 then Liouville Green cannot be used. If it is the case that 
(x — a) 2 V(x) is analytic then one can use the Frobenius method to determine the solutions. This 
procedure is outlined in Q but we do not consider this case in this paper. 



3 Boundary Perturbation 

Theorem 3.1 Let L satisfy the assumptions above and denote the eigenvalues of L by \ n . Now, 
let L e := L |( Q + £ m with Friedrichs boundary conditions and denote its eigenvalues by X n>e . Then 
the following holds: 

A n ,e = A n + c n e p + o(e p ) as e — > 0+ 

where 

p:=2v/c+l/4. 

Proof Let A be an eigenvalue of L. Near x = a we have two linearly independent solutions, 4>i t \ 
and (f>2,\- Since x = a is LCNO both (f>\ \ and <f>2 : \ lie in L 2 near a. We suppose that 6\ \ is the 



recessive solution and so satisfies the boundary condition at x = a - see Corollary 2.2. Also, we 
normalise the functions in order to assume that 

1 ,A(»)~(a:-o) 1 / 2+ V^+^ 

Since A is an eigenvalue </>i,a also satisfies the boundary condition at x — b (either an explicit 
condition if b is R, WR or LCNO or the impicit condition that it must lie in M(c, b) if b is LP). 

Now, we consider the differential equation 

-&(x) + (V(x) - \)£ x (x) = llA (a:). 
At x = a we impose the boundary condition 

[&,6>,A](a)=0 

and at x = b we use the same boundary condition as for the homogeneous equation. 

The corresponding homogeneous equation is solved for the left boundary condition by (f>2,x and 
for the right by </>i,a- These are linearly independent and so we can write down the Green's 
function: 

G(x,y):={ -ft 2 H y \t 1A( } X \ X l V 

I -~<t>iAv)<l>iAx) v> x 

where k := y>i,\,4>2,\](,c) for any c e (a, b). Note that K ^ and is a constant whose value is 
independent of the particular c chosen - see [El p351] for an introduction to Green's functions. 

We write down the function 

h(x) = 0i, aO) + e p c E £\(x) 



where c e is some constant to be determined. We know that this satisfies the right boundary 
condition and next turn to the point x = a + e. We have 

h(a + e) = 4>i t \(a + e) + e p c e £\(a + e) 

and thus choosing 

_ -e~ p 0i,A(a + e) 
£,\(a + e) 

-e~ p ^\{a + e) 

_ « Ia +e <f > 2,\(y)<f ) iA a + £ )4>i,\(.y)dy - j: J a+e <l>i,\(y)<fa,\(a + e)4>i,\{y)dy 
we also have 

h{a + e) = 0. 
Using the asymptotic expansions developed for <pi t \ and (j>2.\ we have the following: 

fa,\{y)4>i,\{y)dy = o(e 2 ) as e —> 0+, 

J a 

b 

2 J-r,(r 2 +P\ 



Thus we have 



4>iAv)4>iAv)dy = \\4>iA\ 2 + °( £ p ) as £ - ° 



«£ P ^ia(a + e) 



k e -p(eV2+p/2 + ( e i/2+p/2) x 

~ H0i,A||i (e 1/2 - p/2 + oCe 1 / 2 -?^)) (i + ( e 2+ P )) 

= tt; — ii9 + °(1) as £ - * + - 
II<pi,a||2 

It is crucial to notice that the leading term in the expansion of c e is independent of the particular 
(pi.x and (j>2. x chosen; i.e. suppose we now took <j)i.\ = A<f>\.\ and 02, a = B<j)2.\, then a calculation 
shows that 

c e = c e 

in an obvious notation. 
Define A e by 

X e := A + c £ e p 

and the function h satisfies 

-h"(x) + (V{x) - Xe)h(x) = -c 2 e e 2 P^ x (x) 



subject to Friedrichs boundary conditions. 

Now suppose A e is not in the spectrum of L e . The operator L £ — X e is thus invertible and we 
have 

\\(L £ -X £ )-^ X \\2= ¥h 



2 F 2p 



( -t 



and so 



ULe ' K) "" II6II2 -\\bh<2&- () 



LetC=||h|| 2 /||6l|2. Then 



c -i = 116,112 _ H6II2 



\h\\ 2 Hi,\ + £ Pc et\h 

< II6JI2 < 2 II6JI2 



1101, aIU -C £ £ p ||^a||2 ||<W,A||2 

for e small enough and thus we have a uniform upper bound for C _1 . 
From (0) we immediately have 

1 C 

> 



dist(A E ,Spec(L £ )) ~ c 2 e e 2 P 

c 2 e 2p 
<^dist(A e ,Spec(L e )) < -^— . 

Thus there exists a point fi € Spec(L e ) such that 

fi = X e + 0{S 2P ) 



2,>\ 



= A + c e e p + 0(e 

WmMh J 

= A + 11 1 K 11 £P + °( £P ) as £ -* ° + • 
II<Pi,a||2 

This conclusion obviously still holds if A e £ Spec(L e ). 

We consider the bottom eigenvalue Ao- The above proves that there is a point /1 £ Spec(L £ ) such 
that 

Ii — Ao < ce p 

for all e small enough. Thus /i eventually becomes smaller that Ai and hence smaller than A„ j£ 
for all n > 1. In other words fi = Xq_ £ eventually. The eigenvalues X n all have multiplicity 1 and 
this argument can be repeated for Ai and so on. This completes the proof. 
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Examples 




4.1 


Example 1 




The theorem applied to the Bessel equation 






-rw + " 2 - 2 1/4 



f(x) = Xf(x) for < x < 1 



< v < 1 



in the case 



yields that fact that 

A„, E = X n + c n e 2u + o{e 2v ) as e -> + . (3) 

This captures the main flavour of the results. In particular it highlights the fact that the rate of 
convergence lies in the interval (0, 2). 

4.2 Example 2 

We now return to the main motivating example. The original problem was two-dimensional. We 
consider the operator 

Hf := -d(x) 2 ^A E 

acting the space L 2 (tt, d(x)~ 2 ~ f ) where J7 is the unit disc in two dimensions. The rotational 
invariance of the domain allows us to separate variables and after some transformations reduce 
the problem to a Schrodinger equation. Here the particular potential is 

V( ^ 2 7-7 2 1 .( a l\ (1 - I ) 27Q 

y (^) = 77^ 7X777 777+ h' 



4(1 - 7 )2 (a -xf V 4; (1_(1_£)«)2 

for x E (0, a). The constants 7 and a lie in the intervals [0, 1/2) and [1, 2) respectively and are 
related by a := 1/(1 — 7). In this case 

V(x) = |^7777-^y2 + (! - 7) 27t V " V4)(« - *r a + 0((a - x)^ +1 )") 
as x — > a— . The results of the paper are then used to show that 

X n ,e = A„ + c„e~ + o{e T ^i) as e — ► + . 
This is in agreement with the abstract result (see 0, Theorem 7.1]). 



4.3 Numerical Example 

We carried out some numerical tests using the program SLEIGN2 - see Jl) for an introduction. 
For simplicity we return to the case of the Bessel equation on (0, 1) and choose v = 0.6. The 
eigenvalues are easily computed using standard computing packages as well as SLEIGN. Next for 
selected values of e we can use SLEIGN to compute the bottom eigenvalue on (e, 1). The results 
are recorded below. 



e 


K 


0.1 


12.6988324 


0.01 


10.8878760 


0.001 


10.7821955 


0.0001 


10.7755528 


0.00001 


10.7751337 


0.000001 


10.7751073 



SLEIGN's quoted tolerances are all smaller than 10 10 . The computed value for Ao was 

A = 10.7751055. 



Assuming that a relation of the form 



A{), ; 



X F = ce d 



making a plot of log(e) against log(Ao. e — A e ) we can estimate the value of 5 from the slope of 
the graph. Performing a least squares fit on the data in the table yields the value 

S = 1.204292602. 

Although this is not a decisive piece of numerical analysis it does provide a reasonable confirma- 
tion of equation (Q) with v = 0.6. 
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